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Abstract 

We develop an effective Gauss-Bonnet extension of Loop Quantum Cosmology, by intro¬ 
ducing holonomy corrections in modified F{Q) theories of gravity. Within the context of our 
formalism, we provide a perturbative expansion in the critical density, a parameter character¬ 
istic of Loop Quantum Gravity theories, and we result in having leading order corrections to 
the classical F{Q) theories of gravity. After extensively discussing the formalism, we present 
a reconstruction method that makes possible to find the Loop Quantum Cosmology corrected 
F{Q) theory that can realize various cosmological scenarios. We exemplify our theoretical 
constructions by using bouncing cosmologies, and we investigate which Loop Quantum Cos¬ 
mology corrected Gauss-Bonnet modified gravities can successfully realize such cosmologies. 
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1 Introduction 


The Big Bang era is one of the less understood periods of the evolution of our Universe, and the 
physics behind this era is still inconceivable. The standard approach that provides a classical phys¬ 
ical description is based on the assumption that our Universe is described by a Friedman-Lemaitre- 
Robertson-Walker (FLRW) metric, which captures most of the present date characteristics of our 
Universe, i.e. its space uniformity, large scale homogeneity etc. However, the classical cosmolo¬ 
gical approach leads inevitably to an initial singularity, which is a rather “embarrassing” feature 
of the classical description, because due to this singularity, the closed time-like geodesics which 
pass from this singularity, have a finite proper length, but no end points to normal space away 
from the singularity. In addition to this, there has been conjectured for some years 111, that naked 
singularities should be well hidden behind horizons, the so-called cosmic censorship conjecture. 
However, no formal proof exists up to date for this cosmic censorship hypothesis, nevertheless it is 
considered a physically appealing hypothesis. From a cosmological initial singularity exists, then 
it belongs to a period of time that quantum physics, or more appropriately, quantum cosmology 
governs the physical phenomena. 

One of the most appealing quantum cosmological theories is holonomy corrected Loop Quantum 
Cosmology (LQC) m, which is promising from various points of view. For an important stream of 
papers and reviews on this fastly developing research topic, consult ||2l, and references therein. In 
the context of LQC, many theoretical inconsistencies or ellipses of the classical cosmological the¬ 
ories, could consistently be explained. Having such a promising theoretical framework at hand, it 
is compelling to investigate theoretical modifications, along the research lines of modified gravity. 
In this context, F{R) extensions of LQC have recently been developed in |[3jlll, where, following 
the idea of l|5l, holonomy corrections are introduced in Einstein frame (EF), because in that frame 
the gravitational part of the Hamiltonian is linear in the scalar curvature and thus the procedure to 
introduce these effects is the same as in General Relativity. These corrections have been applied to 

gravity Q, obtaining a bouncing model free of singularities. On the other hand, inspired in the 
introduction of holomy effects in the theory of cosmological perturbations OTl . an effective way 
to introduce these corrections in F{R) gravity was developed in 13, where the EF formulation is 
not necessary to be used. The main characteristic of it is that it could be applied to other theories, 
where there is not a conformal equivalence with EF, and therefore allowing to introduce these 
corrections to Gauss-Bonnet gravity in an explicit way. Gauss-Bonnet gravity is characterized by 
the Gauss-Bonnet invariant curvature which we denote G in this paper. 

Adopting the research lines of [71, we perform a Gauss-Bonnet an effective extension of EQC. 
Once we have obtained this new theory, we generalize the reconstruction method obtained in [H, 
and we reconstruct some holonomy corrected EQC-F(^) gravity theories. We shall focus our 
reconstruction procedure to realizing mainly bouncing cosmologies. The reason for that special¬ 
ization is that bouncing cosmologies have the very appealing feature of the absence of an initial 
singularity [^ [TOl [TTl [121 HSl [Ml [151 HSl H?! HSl, plus accelerating expansion can consistently 
be described in the context of these theories MMM- These kind of cosmologies are very 
interesting owing to the fact that nowadays they represent the most promising alternative theory to 
the inflationary paradigm, because if at very early times the Universe is nearly matter dominated 
in the contracting phase, the power-spectrum of the modes that leave the Hubble radius during 
this regime will be almost flat ll23 . Moreover, is has been shown in ll2^ that bouncing cosmolo- 
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gies provide theoretical values of the cosmological parameters, such as the spectral index and its 
running, that fit well with recent Planck data Il25l . 

But what is the motivation to introduce holonomy corrections in an F(^) gravity in the first 
place. We need therefore to discuss the motivation to include LQC effects in F{Q), since Gauss- 
Bonnet as a kind of effective limit of string theory, and the Einstein-Hilbert LQC corrected gravity 
is finite, so why introducing LQC effected F{Q) gravity. The answer is that, since the non-LQC 
F(0)-gravity contains finite time singularities, which in some cases can be crushing type singular¬ 
ities. So a necessary extension of standard F{Q) gravity to include LQC effects is necessary to see 
if any new effects that the LQC corrections might bring along. Lor example, recently in Ref. |[26]l . 
it has been pointed out that for a Type IV singular bounce, near the bounce point, in which case 
all the modes that are relevant for today observations are well inside the Hubble radius, the power 
spectrum of primordial curvature perturbations is not scale invariant. Thereby, a quite important 
task is to investigate whether the spectrum of primordial curvature perturbations becomes scale 
invariant by including LQC holonomy corrections in the classical F{Q) gravity. The same case 
applies for the Einstein frame LQC corrected F{R) gravity, with a thorough investigation of this 
case being performed in @. Hence, the LQC corrected F{Q) is necessary for the aforementioned 
reasons, and also in order to see whether the singularity structure of the cosmological dynamical 
system is changed at some extent. 

The paper is organized as follows: In section 2, using the method of Lagrange multipliers, we 
construct the classical dynamical equations in Gauss-Bonnet gravity, for a flat LLRW geometry. In 
section 3, we review the way to introduce holonomy corrections in standard Einstein-Hilbert gen¬ 
eral relativity, for the particular case of the flat LLRW geometry. In section 4, we extend LQC to 
Gauss-Bonnet gravity obtaining a holonomy corrected Lriedmann equation, which contains all the 
dynamical information of the system. Section 5 is devoted to explain the reconstruction method 
in holonomy corrected Gauss-Bonnet gravity. Sections 6 and 7 contain examples of reconstruc¬ 
tion, for power law, bouncing exponential models and generic bouncing models respectively. The 
conclusions follow at the end of the paper. 

2 Gauss-Bonnet F{Q) gravity ”a la Ostrogradsky” 

When one considers the flat LLRW geometry, the Lagrangian of Gauss-Bonnet F{Q) gravity in 
the vacuum is given by (in units h = c = SttG = 1), 


C{V,V,V) = ^ViR + F{g)) 


( 1 ) 


where V = is the volume and g denotes the Gauss-Bonnet curvature, which can be expressed 
in terms of the volume and its higher derivatives. 



( 2 ) 


Ostrogradsky’s idea to obtain the Hamiltonian from a Lagrangian containing higher order derivat¬ 
ives, is epitomized in the introduction of a Lagrange multiplier term in the Lagrangian, which we 
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denote n, in the following way |[27l |71 : 


c,{v,v,v,g) = -v{R + F{g)) + 


8^2 / V 


2V 


>2 


3^2 I 3^ 9V2 


-g 


(3) 


Maximizing with respect to the Gauss-Bonnet invariant g yields ^ = ^VF'{g). In order for 
the second derivative of V to be removed, we subtract from the Lagrangian the following total 
derivative term, 


(4) 


which does not change the dynamics of the system, and by replacing the Lagrange multiplier /r by 
its value iVF'(g), we finally obtain. 


civ,v,g,g) = -^ + ^VFig) --VF\g)g - ^F'\g)g 


(5) 


As it is obvious by Eq. ([5]), we can see that the Lagrangian ([5]) depends on the variables {V, g) and 
their first derivatives. We can therefore obtain the corresponding canonically conjugate momenta, 
which are. 


Pv 


_dc{v,v,g,g) 2v 4v^ 


dv 


3E 9^2 


F"{<g)g, 


pg = 


dc{v,v,g,g) 4 i /3 


F"{g), 


dg 27 E 2 

and consequently the classical gravitational part of the Hamiltonian becomes, 


( 6 ) 


F- 

2/3 


graviv,g,pv,pg) = Vpv + gpg - c{v,v,g,g) 


= V (^ (aF'(g) - F(g)) 


V [4F"(g) J 


\4F"{g) J 


(7) 


The Hamiltonian constraint T-LgraviF, g^Pv^Pg) = 0 leads to the well-known modified Friedmann 
equation in the vacuum for non-LQC corrected F{g) gravity, which is. 


6772 ^ 24H^gF"{g) - gF'{g) + Fig) = o. 


( 8 ) 


Note that in ordinary Einstein-Hilbert gravity, where Fig) = 0, the canonically conjugate mo¬ 
menta are equal to. 






(9) 


and therefore, the classical gravitational part of the Hamiltonian takes the following simplified 
form. 


n 


graviV,Pv) = -^pIV = - 377 ^. 


( 10 ) 
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Remark 2.1. Note that, the Hamilton equations V = and Q = are simple iden- 

titles. The equation pg = - is equivalent to the modified Friedmann equation (|S|), and 

dl~L 

the equation p'y =- which corresponds to the modified Raychauduri equation in F{Q) 

gravity, could be obtained taking the derivative with respect to the cosmic time of the modified 
Friedmann equation. a consequence, in the vacuum the dynamics in F{Q) gravity is only mod¬ 
elled by the modified Friedmann equation (jS]), and when one considers matter one has to include 
the conservation equation p = —3H{p + P), being P the pressure. 

Before ending this section, we need to stress an important issue: Since in principle there 
are infinitely many canonical transformations, this means that the modified F{Q) gravify could 
be formulaled using infinilely many sefs of variables (fwo coordinates and fheir corresponding 
canonically conjugate momenfa). Nofe fhaf some of fhese sefs of variables will be meaningless 
physically speaking, because fhey are builf using a combination of bofh coordinafes and momenfa, 
giving new quantifies wifh a very difficulf physical inferprefafion. Moreover, since fhe infroducfion 
of holonomy effecls critically depend on fhe sef of variables used, in effecl fhere are infinitely 
many ways fo infroduce holonomy effecfs in modified F{Q) gravify. Consequenfly, fhere will be 
infinifely many differenf effective holonomy corrected Friedmann equations in F{G) gravify. 


3 Introduction of holonomy corrections 


In fhis secfion and in order for fhe arficle fo be mainfained self-complefe, we shall describe fhe 
technique of infroducing holonomy corrections in sfandard Einsfein-Hilberf gravify. Assuming a 
flaf FLRW geomefry, firsf of all one can consider fhe variable = —^Pv = 7^ ( 1291 ). where 
7 is fhe Barbero-Immirzi parameter. In terms of fhe paramefer fhe Hamilfonian (fTOt becomes 
T~Lgrav{y, fi) = However, in LQC, due fo fhe discrete nafure of space, fhe quanfum 

operafor /3 is nol well defined (see for insfance l30l or ISTl for a sfafus reporf on fhis issue). Then, 
in order fo build fhe quanfum fheory, fhe gravifafional part of fhe Hamilfonian musf be redefined. 
To be precise, we will consider holonomies of fhe form hj{\) = where Oj denote fhe 

Pauli mafrices and A is fhe square roof of fhe minimum eigenvalue of fhe area operafor in loop 
quanfum gravify. Since /3^ does nof have a well-defined quanfum operafor, in order for a consisfenf 
quanfum Hamilfonian operafor fo be consfrucfed, an almosf periodic funcfion fhaf approaches 
for small values of is needed. This can be done using fhe general formulae of loop quanfum 
gravify fo obfain fhe holonomy corrected Hamilfonian, which is equal fo. 


nhoi,grav{V,/3) = [hi{X)hjiX)hr\X) 

xh-\X)hk{X){hl\XfiV} . 


( 11 ) 


The Hamiltonian (fTTl) captures the underlying loop quantum dynamics, and for a detailed account 
on this issue see for instance 132]]). 

A simple calculation ll^ l34l [35l shows that the Hamiltonian of Eq. (fTTl) acquires the simple 
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form 


nhol,grav{V, /3) = (12) 

which indicates that effectively, holonomy effects can be introduced by explicitly performing the 
replacement /3 sm(\i3) equivalently pv —)• 

In order to obtain the holonomy corrected Friedmann equation, in principle it is necessary to 
use the full Hamiltonian, i.e., the Hamiltonian that contains both, gravitational and matter parts, 

-HhoiiV, /3) = + pV, (13) 

Then, by combining the Hamilton equation, 

1 dUhoiiV, P) ^-,^ sin(2A/3) 

2 9/3 2 A 7 ’ ^ ’ 

with the Hamiltonian constraint T-LhoiiV, P) = 0, we obtain the well-known modified Friedmann 
equation Il36]l . 


= 


P 

3 



(15) 


where pc = is the so-called critical density. Practically this density is what measures the 
strength of the loop quantum effects in cosmology, so if it is infinite, and for finite energy densities, 
the holonomy corrected Friedmann equation will be reduced to its well known form in standard 
cosmology. 


4 Holonomy corrected F{Q) gravity 

Recently, holonomy correction are introduced in F{R) gravity using the property that F{R) grav¬ 
ity in the Jordan Frame (JF) is equivalent to General Relativity in the Einstein Frame (EF). Then, 
introducing, as in standard EQC, holonomy corrections in EF and comming back to the JF one 
obtains what is named as Eoop Quantum F{R) gravity 131IH. 

Unfortunately, Gauss-Bonnet gravity is not equivalent to General Relativity in any frame. 
Then, to introduce holonomy corrections, other strategy must be used. What we will set up is 
an effective EQC theory for F{Q) gravity based in this key point: As we have seen in previous 
Section, stardard holohomy corrected EQC, could essentially be obtained from the replacement 
j3 Qj. equivalently pv —>• Therefore, when there is not an established way to 

perform a kinematical loop quantization of the phase space, as when one deals with cosmological 
perturbations, in order to incorporate loop corrections, an effective Hamiltonian is built adopting as 
a prescription the replacement of the Ashtekar connection by a suitable function of that connection 
IITTI . Then, with the same spirit, we will introduce holonomy corrections in general F{Q) using 
the following recipe: 
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To introduce the holonomy correction in general F{Q) we will adopt the following recipe: In 
analogy with the standard Einstein-Hilbert gravity case, where F{Q) = 0, we will replace the 
momentum that in Einstein gravity case corresponds to —2^ by por example, if the 

variables {V,G,pv,Pg) are used, we make the replacement pv —)• in the Hamiltonian 

of Eq. (|7]l. It is conceivable that this way of introducing holonomy corrections, critically depends 
on the set of variables used to formulate the F{Q) theory, which means that the use of other 
canonically conjugate variables will lead to different corrections. In the case that matter fluids 
are taken into account, considering the Hamiltonian (|7]l and its corresponding matter part, namely 
'Hmatter = pV, upou making the replacement, pv —>• — we finally obtain the following 

Hamiltonian, 


6sin(A/3) 


T-ihoi{V, G,pv,Pg) 


H \4F"ig) J 


( V^pg 
\4F"{g) 


1/3 Y 

+ -{gF'ig)-F{g))+pv. 


(16) 


The corresponding Hamilton equations are equal to. 


7 dUhoi ^ _ dTihoi 
2 5/3’ dpg 


(17) 


which together with the Hamiltonian constraint T-LhoiiV, G,Pv,Pg) = 0, lead to the following 
equations: 


g = 


1 


2F"(g) 


H = — cos(A/3)p| 

At j 


Pg +Pg 


^-4 sin 


3~| + ~l ^ 1 ^ ^ ^ 


(18) 


where for notational simplicity, we introduced the variable pg = 

Before dealing with these equations, we will check that in standard Einstein-Hilbert gravity 
where F{g) = 0, we can obtain the holonomy corrected Eriedmann equation of Eq. (fTSl) . In 
Einstein-Hilbert gravity, the equations (fTSl) take the following form. 


Pg +Pg 


H = — cos(A/3)p| 
,-|sin(A/3) 


3p| -h 


3 6sin(A;0) 


At 


At ° 

+ p = 0. 


(19) 


2 

The second and third equations lead to the relation = f, and the first and second equation can 
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be written as follows, 


sin^(A/3) = 1- 2 " = 1 “ 

Pg 


P 


sin^(A/3) 



P_ 

Pc 


( 20 ) 


which eventually lead to the Friedmann equation (fTSl) in LQC. Coming back to the general equa¬ 
tion (fT^ . these thr'ee equations have to be used to obtain a general relation of the form G{H, Q,Q, p) 
0 which will correspond to the modified Friedmann equation in Gauss-Bonnet F{Q) gravity. This 
equation in conjunction with the energy-momentum conservation equation p = —3H{p -|- P) and 
the equation of state P = P{p), we the equations that will depict the dynamics of the system. 

The combination of the first and third equations of (IT^ leads to the relation, 


12p,(p| - H^) = Q {QF'{Q) - F{Q)) + p + 3p|) , 
2 

which allows us to isolate as a function of H, Q and p, giving as a result 

4 Ap,-A{g) (^ 1^ A^G) + 46f^ \ 

6 V (4p,-A(a))2 


( 21 ) 


( 22 ) 


where the notation A{G) = QF'[Q) — F{Q) + 2p has been introduced. 

On the other hand, the combination of the three equations in (fT^ leads to the dynamical 
equation that corresponds to the modified Friedmann equation in LQC-F(^) gravity. 


= f 5 j - ^(^)) - f • (23) 

To sum up, the dynamical holonomy corrected equations in Gauss-Bonnet gravity are given by the 
modified Friedmann equation in LQC-F(^) gravity and the conservation equation. 


= f ^ - no)) - f 

p = -3H{p + P), (24) 

2 

where pg also appears in Eq. (1221) and pc is the critical density introduced at the end of the previous 
section. 

Here, as in standard F{G) gravity, only two equations are needed to define the dynamics. All 
the Hamiltonian equations obtained from (fT^ are equivalent to (1241) . 

Note finally that, in order to solve Eq. (1241) . for a barotropic fluid an equation of state of the 
form P = P{p) is needed. In contrast, when one deals with a canonical scalar field tp, one has 
p = -h V(p) and P = - V(p). 






5 Analysis of the theory F{Q) with LQC 

In this way we shall investigate in which way we can have some LQC corrections in modified 
F{Q) gravity, at least the leading order LQC corrections. These corrections would materialize the 
first deviations from the classical non-LQC F{Q) gravity. In order to find these, we rewrite Eq. 
d2^ in the following form, 

2^{F"{Q)Q?Pg{Q) - 2pc + A{g) = 0, (25) 

where the pai'ameter pg{G) stands for, 

Pg{Q) = ^ (4pe - A{Q) - - A{Q) - . 

Note that when pc —)• oo, the quantity pg{G) can be represented as a series 

P.(e) = + T (Am +i + ^ (.4(8) + ^■ 

Then, Eq. (1251) can be written as follows, 

576H\F''{g)gf - {A{g) - + 

+ e ^ + 1152H^{F''{g)g?) + 

^, {Aig) + 6H^)\ 

2304iT4 

X (^A^{g) - 6A^{g)H^ + 432/7® - 576H^{A{g + 30H‘^)){F"{g)gf^ + ... = 0. (26) 

It is easy to see that in the limit pc ^ oo this equation reduced to the Eriedman equation for F{g) 
theory (see 1381). Therefore, we define the parameter e = 1/pc, which as pc —)• oo, it takes small 
values, and we seek a solution of this equation by performing a perturbative expansion in terms of 
this parameter. 

In order to provide a consistent solution to the Cauchy problem or a boundary value prob¬ 
lem, for the differential equation (1251) . we shall perform a perturbative expansion in terms of the 
parameter e, as e —)■ 0. This expansion takes the form, 

CX) 

F{g) = Y,FFk{g). (27) 

k={) 

Then, the differential equations for each order of the expansion are given below, 

2 m^F”{g)g - gFl^{g) + Fo(g) + 6/7^ - 2 p{t) = o, ( 28 ) 

2ah^f”{Q)G - GF[{g) + Fi{g) = + 2HF”{g)gf{i + QHF”{g)g), m 
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2AH^F”{G)G - GF^iG) + F2{G) = 

= -9H\3H{1 + 2HFi;{G)Gf{A + 39FF^'(0)^(1 + 2HFi;{G)G))+ 

+ 40(1 + 2HF”{G)G){f> + imF”{G)G)F”{G)). (30) 

where we included only the first three orders of the perturbative expansion of the solution. An 
important comment is in order: In principle, it is expected that the entire set of dynamical variables 
should be expanded in terms of the parameter e, however if this is done for the Hubble rate, this 
would lead to a very restricted final form of the solutions. What is considered as a perturbative 
solution in the case at hand, is the full form of the F{G), and the Hubble rate will be considered 
that it has a specific form. So we actually provide a perturbative method of reconstruction of 
the F{G) gravity that generates such an evolution. In fact, the parameter pc, and effectively e, 
is justified only by the LQC context, and it is not expected to appear in the Hubble rate or in 
the matter energy density p{t). The latter two physical quantities will be determined from the 
beginning, and in the later sections we shall investigate which F{G) LQC-corrected gravity can 
generate the cosmological evolution corresponding to the given Hubble rate and matter energy 
density. Therefore, e enters the perturbative expansion only via Eq. dTTl) . If it is possible to 
construct a solution of equation (|2^ for the function Fq{G), the remaining terms of the expansion 
Fk{G), k > I are determined by solving linear differential equations. The initial or boundary 
conditions for the functions Fk{G) can be obtained by taking into account the expansion (IZ/I) and 
the conditions for the equation (1251) . 

Note that from the equations defining the function Fk{G), it follows that the corrections to the 
Fq{G) are not related to the distribution of matter in the universe. In the following sections, we 
shall exemplify our results by using illustrative examples, emphasizing in bouncing cosmologies. 
We shall assume that all matter fluids are absent i.e. p{t) = 0 , and therefore we are interested in 
vacuum modified LQC-corrected F{G) theories. 

6 Power-law Cosmology from LQC-corrected F{Q) Gravity 

We start of our analysis, by studying a power-law cosmology in the context of LQC-corrected 
F{G) gravity. Consider an evolution of the Universe for which the scale factor is of the following 
form, 

a{t) = a e M, n G N. 

Lor this model the Hubble rate reads. 

Whence, it follows immediately that, 

Hc=- ^ _ 

2V3V2re - I v/3v^V2n - I 

We assume that the initial conditions for the Cauchy problem of the differential equation (l25l) . are, 

F(0o)=7i, F'{Q^)=12-, 
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where 71, 72 are constants, 0 < < +00. Then, for the system of differential equations for the 

functions Fk{Q), we obtain the following initial conditions: 

Fo{Go) = 7 i) -^o(^o) = 72, 

Fk{go) = 0, F^(eo) = 0, A: = 1,2,.... (31) 

The general solution of the differential equation for the zero order unperturbed function Fq{Q), 
namely Eq. (1281) . is of the following form: 


FoiG) = CiG + 


2V3V^y/n(2n - 1) 

2n + 1 


where ci and C2 are arbitrary integration constants. 

In order to further simplify our analysis, we can choose 71 and 72 in such a way, so that 
Cl = C2 = 0, without great loss of generality. In effect, Fo{G) reads. 


^ 2V3Vgy^n{2n-l) 

I 1 


(32) 


Having Fo{G) at hand, we can easily proceed in finding the first order correction of the perturbative 
expansion, namely the function Fi (G), which is of first order in the e expansion. Considering Eqs. 
(OTI) and (l32l) . the Cauchy problem for (|2^ is the following: 


^^F^iG)-GF[iG) + F,iG) 


12Gn^in-l) 

(1 - 2n)(l + 2n)3’ 


Fi{Go) = 0, 


= 0 . 


Therefore, one analytic solution of the Cauchy problem is: 

F.(e) = .36 + + (r Itntf 2n)^ ■ 

12(1 — n)n^\n.Go 48(1 — 

~ (l + 2n)3(3 + 2n)’ “ (1 + 2n)3(3 + 2n)2 

Note that when n = 1, the Cauchy problem has the trivial solution, Fi{G) = 0. Iteratively, we 
can calculate the second order correction in the e-expansion. This can be easily done, since we 
have the explicit form of the function Fi{G) at hand. Combining Eqs. (OTI) and (021) and also by 
taking into account the explicit form of the function Fi{G), the Cauchy problem for the differential 
equation (OOl) . takes the form : 


-^-F!^{G)-GF!,{G)+F2{G) = 

^ V^n3/2(2n + 3)(5n - 2) = 3 _ 2 nV 4 3V3n^/^{32n^ - JSv? + 51n - 20) 3/2 

“ 4v^2^)r^(2n + 1)2 (2n- l)3/2(2n + l)5(2n + 3) 

-^2(^0) = 0, ^2(^0) = 0. 
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It is easy to show that the solution to this problem has the following form, 


F,(g) = c,g + + ,, + 12n^ - 23n + 6) 


+ 


+ 


6V3n9/2(32n3 - 78n^ + 51n - 20) 3/2 
V2n-l(2n + l)5(2n + 5) 


where C 5 , cq are constants, which is easy to specify from the initial conditions given above. In 
principle, by continuing this iterative process, we can find all the higher order corrections of the 
perturbative e-expansion. 

In conclusion, for the case n = 1, the resulting approximate solution of Eq. (l25l) . has the 
following form. 




2VU 2(2VG^ 1 

^ 45 ^ 3 ^ 315^3 01 /^ 



-I- o(e^). 


(33) 


where we took into account only the first two corrections of the perturbative e-expansion. As 
a final fask, we invesfigafe the discrepancy between the approximate solution (l3^ and the one 
appearing in Eq. (1251) . Eor this purpose we substitute (1^ in equation (1251) for this model and as 
a result we obtain the function g{G,s) = o(e^), e —)• 0. We have to note that the approximate 
solution is defined in fhe neighborhood of poinf G = Gq. 


1 Bounce Cosmology from LQC F{Q) Gravity 

One appealing cosmological scenario which is alfernafive fo fhe sfandard inflafionary description 
is fhe bounce cosmology scenario ifTOl fTTl [T2j [T3l [141 [El [El [El [El- According fo fhis, fhe 
Universe confracfs until a minimal radius is reached, where if bounces off and sfarfs fo expand. In 
the bouncing cosmology context, the initial singularity is avoided and this feature is mainly what 
renders the bouncing cosmologies so appealing, in comparison to the inflationary picture. Also, 
it is possible to consistently describe early-time acceleration |[20ll^l22l . so in conjunction with 
the singularity avoidance, we have a very appealing candidate theory for our Universe’s evolution. 
In the context of modified gravify, bouncing cosmologies are easily realized |[T71[T8l[T9l, wifhouf 
the peculiar requirements that the standard Einstein-Hilbert imposes in order for a bounce to occur 
ifTOil . In this section we shall investigate which EQC-corrected F{G) gravities can realize bouncing 
cosmologies, using well known bouncing cosmologies paradigms. Special emphasis shall be given 
on the exact form of the EQC-corrected F{G) which describes the bounce near the bouncing point. 
We address the problem of finding the EQC-corrected F{G), using the perturbative expansion we 
used earlier. 

7.1 Exponential Bouncing Models 

We start off our analysis by studying some bouncing cosmologies with exponential scale factor. 
Before we get into the details of each model, it is worth recalling in brief the conditions that define 
a bounce cosmology. Eor defailed accounfs on these issues, see for example |[T0l[m[T^[T^[T4l[T5]l . 
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As we already mentioned, a cosmological bounce is described by two evolutionary eras, a 
contraction and a subsequent expansion. In the contracting era, the Universe’s scale factor a{t) 
decreases up to a point, say at t = at which a minimal radius is reached. This means that during 
the contraction, the derivative of the scale factor is negative d < 0, and at the minimal radius point 
the derivative of the scale factor is zero, that is d = 0. At this point, the Universe bounces off and 
starts to expand, so that d > 0, until a finite time singularity is probably reached. We shall study 
two bouncing models, the scale factor of which contains exponential functions of the cosmological 
time t. 

7.1.1 Symmetric Bounce Model 

In this section we study a symmetric bounce model, the scale factor of which is equal to, 

a{t) = , (34) 

where a is a positive arbitrary parameter. The scale factor (1341) describes a cosmological bounce, 
for which the bouncing point is at f = 0. Indeed, this can be verified in Fig. [T] where we ploffed fhe 
fime-dependence of fhe derivative of fhe scale factor d. As if can be seen in Fig. [T] fhe function d is 
negative before fhe bouncing poinf f = 0, equal fo zero af fhe bouncing poinf and posifive affer fhe 
bouncing poinf. Whaf inferesfs us fhe mosf in fhis section, is fo invesfigafe which LQC-correcfed 



t (sec) 


Figure 1: The derivafive of fhe scale factor a{t) as a function of fhe cosmic lime t, for fhe bouncing 
cosmology a{t) = e“* , wifh a = 1. 

F{G) fheory can generate fhe bounce (l34l) . near fhe bouncing poinf t = 0. If is conceivable fhaf 
whaf we need to examine is fhe limif t —)• 0 of fhe differenlial equations (1281) and (l29l ) presenled 
in fhe previous section. We shall be interested in finding fhe firsf fwo functions of fhe perlurbalive 
e-expansion of fhe previous section, namely Fq{Q) and Fi{Q), assuming absence of any mailer 
fluids. The Gauss-Bonnel invarianl for fhe scale facfor (l34l) reads, 

G = (l + 2t^a) , (35) 

and in fhe limif t —)■ 0, fhis is approximately equal to, 

g ~ . (36) 
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Therefore, Eq. (l36l) can be explicitly solved with respect to the cosmic time t, 


8V3a3/2 


( 37 ) 


Notice that the small t limit is equivalent to the small Q limit, so we shall take this observation into 
account in the rest of this section. In view of the above approximations, the differential equation 
that yields the function Fq{Q), namely Eq. (1281) . at the vicinity of the bouncing point f = 0 is 
equal to. 




(38) 


which can be solved analytically, with the solution being of the following form. 


Fo{g) = gci 


y/gy/^^C 2 ln[^] + 2^/gy/A&a^ In \/48 q;2 


24a2 


8a 


(39) 


where we kept only the leading order terms in the limit ^ 0. Having Fo(0) at hand, we can 

easily obtain the function Fi{g), by solving the second differential equation given in Eq. (1291) . 
In order to find an analytic approximation of the solution Fi(^), we need to find an approximate 
form of the term, 

c = - 1877^(1 + 2 HFi;{g)gf{i + QHFi;{g)g ), ( 40 ) 

appearing in the Fi{g) differential equation, where for notational simplicity we denoted the term 
with C. The resulting expression of C is quite lengthy and can be found in the Appendix, so here 
we quote only the simplified expression for C, where we will keep only the dominant terms in the 
limit 0 —)• 0. In this limit, the term C reads. 


C = AiVg, 


(41) 


where Ai stands for. 


-Ai = - 

9 


+ — 


|v^ln 


4v/3v^l 9V3C; 


2 Va^ In 


4\/3\/^ 


128v/3a 


64a2 


32a 


(42) 


H-v/Sv/^ln 4\/3\/^ 

16 L 


■1 3 


Consequently, the differential equation that yields the solution for Fi{g), namely Eq. (I29l ). be¬ 
comes. 


1 


a^ 


48 + -^) F”{g) - gF[{g) + Fi{g) - AiVg = o, 


(43) 


which can analytically be solved to yield. 


Fi(g) = e-S- gCi + — 


(44) 
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In conclusion, the final form of the small e limit of the F{Q), which generates the bounce (1^ 
near the bouncing point, is approximately equal to. 


F(g)~Fo(g)+eFi(g)~ 


QCi- 


\/^\/4^C2 -Gln[g]-\-2^/gV4^ln 


24a2 


+ 


e e 


-gci + 


8a 

12 a^ 12 y^ 

e e e gC2' 


(45) 


Recall that the variable e is equal to l/pc, so practically the limit e —)• 0 measures the differences 
of the resulting LQC-corrected F{g), in the limit where the LQC effects are strong, and therefore 
the LQC effects are contained in the term which is linear to e. Indeed, the bounce is expected to 
occur when p = pc, and at that point, ^ oo, and also, since pc is quite large, the parameter e is 
quite small, and in the expansion, the most dominant term of the perturbative expansion, is of the 
order ~ e. 


7.1.2 Hyperbolic Cosine Bounce Model 

Another exponential bouncing model is described by the following scale factor, 

a{t) = cosh(Af), A > 0 
For this model the Hubble rate is equal to, 

H{t) = Atanh(Af), g{t) = 24A^ tanh^(Af). 


Whence it follows immediately that. 


= 


24A4’ 


HQ 


2\^g 



Consider the following initial conditions for the Cauchy problem of Eq. (1251) . 

FiGo) = 7i, F'{go) = 72, 

where 71, 72 are constants, and also ^0 satisfies 0 < < 24A'^. Then, for the functions Fk{g), 

the following boundary conditions hold true: 

-^o(^o) = 7i) -^o(^o) = 72, 

Fk{Go) = 0, Efc(^o) = 0, A: = 1 , 2 ,.... (46) 

The general solution of the unperturbed differential equation (|2^ has the form: 

Fo{g) = ciG + C 2 -\/ 0 \/ 24 A'^ — G— 
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where ci, C2 are constants of integration. Note that this solution is defined for 0 <G < 24A^. It is 
not difficult to find the values of the constants ci and C2: 


Go-l 2 X^ , 2 AX^-Go , I , f 

Cl = 7i—rrmT;-1" 72—-1- In ' - 


l2X^Go 


12A4 ^2A2 \ 2 Vexy' 


V24A4 - Go V24A4 -GoVG^ 1 , 

= —m3- 


^ 2\/6A^ + y/Go 
I V24A4 - Go I 


V 


Note that in the limit Go 24A‘^: 


7i 71 

C2^-^. 


One can consider the first-order correction of the perturbative e-expansion, namely Fi{G)- By 
considering Eqs. (I46b and (1471) . the Cauchy problem for the differential equation (l29l ) becomes: 

2^2 (1 - FUG) - GFiiG) + F^{G) = 

- 2^XUfUG) - 6X^)UG{G - 24X^)FUG) - 2A2), (48) 

F,{Go) = 0, Fi{Go) = 0. (49) 

It is easy to show that the solutions of the homogeneous differential equation 

yiiG)=G, y 2 {G) = VGV 24 X^-G. 

Then the solution of the inhomogeneous equation (|4^ can be written as: 


are: 


Fi{G) = C 3 G + c^'JG'J 24A'^ — GF 

+ 23^ /" 24A4)F''(g) - 2X^){G{G - 24XUfUQ) - ^X^fdG- 

^23l%o~ ^ - 24A")F^^(g) - 2A^)(g(g - 24A^)F^^(g) - 6A^)^dg, 


where C3, C4 are constants of integration. We can explicitly calculate the integral by using a new 
variable, 


s = arctan 


y/G + 2V6A^ 
V24A4 - G 


TT TT 

— < S < —. 

4 “ 2 
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Then, the first correction Fi{Q) can be written as follows, 


Fi{s) = 12A'^(2 c 3 cos^ 2s — C4 sin4s)+ 

+ (1 + 4s^ + 76(s + C2 A^)^ + 72c 2A'^ ln(24A^)) cos^ 2s— 

9 9 

— 7C2A® (9 In 2 + 11 ln(24A^)) sin4s — -A^ arctan(cot 2s)(8c2A^ cos^ 2s + sin4s)+ 

4 8 

+ ISA"^ (4(7,A'(s) + 2,S'(s)) cos^ 2s + (18,A(s) - 18,S(s) + 5,A"(s)) sin 4s) - 
— 18A^(s + C 2 A^) ((18Al'(s) — 18;S'(s) + 5 ln(sin 2s)) sin 4s + 4(7Al''(s) + 2;S"(s)) cos^ 2s) + 
+ 18A^(s + C 2 A^)^ (9(Al"(s) — B"{s)) sin 4s + 2(2 In | cos 2s| + 7 In sin 2s) cos^ 2s) — 

— 54A'^(s + C2A^)^ (cos 6s esc 2s + sin4s In | tan 2s|), (50) 


where the functions Al(s) and B{s) are defined as follows: 


,A(s) 


11 q 1 O ^ 

-,“ + -8=>ln2a + ^(-l)’ 


n=l 


o4ti—1 

^ _p ^ 2 fc +3 

k{2k + 3)]^ 


n=l 


22”-4(1 - 2^^) 
fe(2A: + 3)! 


2s)2fc+3^ 


and B 2 k are the Bernoulli numbers. Note that 


Al'''(s) = lnsin2s, B'''{s) = ln|cos2s|. 

From the initial conditions (|4^ . it is easy to find the coefficients C3 and C4, but due to lengthy 
and complicated expressions, we do not present the final solution to the Cauchy problem. 

Before closing this section, we need to note that the following limiting cases hold true, 

lim Fq{Q) = 0, lim Fq{Q) = 24ciA4. 
g^o+ g^2A\^- 


lim Fi(^)=0, lim Fi{Q) 
P^>0+ P^24A4- 


|A4 (1 + 187r2 ln(24A4) + 224,4' (f)) , C 2 = 
+00, C2 ^ 


This implies that for the case C2 / fhe first order correction Fi{Q) is defined only in a 

neighborhood of the point Q = Qq. However, the case C2 = — ^ should be considered a separate 
way. It can be obtained, for example, in the formulation of the boundary value problem of Eq. 

dB]): 

lim FiQ) = 0, lim F(Q) = 0. 

g^0+ 5^24X4- 


At the same time it is easy to show that. 


Ci = 0, C2 = -^, C3 = -^(l + 187r2ln(24A4) + 224Al'(|)), 


C4 G M. 
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Note also that in this case, the first order correction Fi{Q) can take finite values on the whole 
interval, where the function Fq{Q) is defined, which is the following interval, 

< C, 0 < a < 24A^ 

where (7 is a finite constant. 

Similarly, one can build corrections any other order for the model a{t) = cosh(Af), but the 
expressions can be more complicated so we confine ourselves to the first two orders. 


Fi{Q) 

Fo{G) 


7.2 Power-law Bouncing Model 

As a final example, we shall study another bouncing cosmology which is described by the follow¬ 
ing scale factor, 

a{t) = {t-tsT, (51) 

which is related to certain ekpyrotic models ifTTlfTSl . The parameter a is a real positive number 
which for the purposes of this section we choose it to satisfy 1 < a < 5, for reasons to be clear 
later on in this section. The cosmological evolution described by the scale factor (ISTT) perfectly 
describes a bounce, meaning that before the bouncing point, we have d < 0, after the bouncing 
point d > 0 and at the bounce d = 0, as can easily be checked by looking at Fig. |2j where 
we plotted the time dependence of the function d, for a = 4/30 and tg = 10 sec. Since for 



-40 -20 0 20 40 

t (sec) 


Figure 2: The derivative of the scale factor d{t) as a function of the cosmic time t, for the bouncing 
cosmology a{t) = {t — with a = 4/3 and tg = 10“^^sec. 


general values of a, it is quite difficult to analytically solve the differential equations (l29l ). as in the 
symmetric bounce case, we shall investigate here which LQC-corrected F{Q) gravity can describe 
the bounce (|5T1) near the bouncing point t = tg, which can be arbitrarily chosen. Notice that when 
t ~ tg, the expression x = t — tg tends to zero, and this observation shall be useful in the following 
analysis. The Gauss-Bonnet invariant Q, for the scale factor (l5Tl) reads. 


24(-l + a)a3 


(52) 


'Notice that in principle if a is not appropriately chosen, the scale factor might turn complex, so extra attention 
should he payed on this issue. 
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so as t tg, the variable Q increases. Therefore the limit x —)• 0 corresponds to the limit ^ oo, 
and we shall use this correspondence in the analysis that follows. Solving Eq. (1521) . with respect 
to t — ts, we obtain, 


t-ts 


23 / 431/4 (_«3 +^ 4 ) 1/4 
gl/4 


(53) 


and since 1 < a < 5, no inconsistency related to complex cosmological times occurs. The 
differential equation that yields the Fq{Q) gravity is therefore equal to, 


- gF'(g) + FoiG) + vVG = 0, (54) 


where we have set for simplicity V and Bi to be equal to. 


V = 




Bi = - 


■\/ (—1 + a)a^ a — 1 

The solution to the differential equation (1541) is equal to, 

, , iV^/G 

^O(0} = Cl + GC 2 . 

— I + Ol 


(55) 


(56) 


Finally, we calculate the first order correction of the LQC-corrected F{G) gravity, namely the 
function Fi{G), which easily follows if we use the analytic form of Fq{G) given in Eq. (1561) . 
However, since the exact solution is quite complicated and lengthy, we can approximate the res¬ 
ulting expression by recalling that the limit x —)• 0, corresponds to ^ 00 , so the larger term 

dominates in the expressions. For 1 < a < 5, the most dominant term is of the order ~ G, so 
F{G) is approximately equal to Fo(^) ~ C 2 G, and the resulting differential equation that yields 
the solution for Fi{G), namely Eq. (I29l) . becomes. 


BiG^F”{G) - GF[{G) + Fi{G) - 

which can analytically be solved to yield, 

Fi{G) = G^^Ci + GC2 + 


?,G(F 


= 0 , 


4 (—a3 + a^) 

3Ga{-Bi-ln[G]+Biln[G]) 


4(-l+Hi)2(-l + a) 


(57) 


(58) 


8 Conclusions 

In this paper, we extended the holonomy corrections formalism of EQC to Gauss-Bonnet F{G) 
modified gravify fheories. Specifically, upon using fhe mefhod of Eagrange mulfipliers, we con- 
sfrucfed fhe classical dynamical cosmological equafions in Gauss-Bonnel gravify, for a flaf FERW 
geomefry. In addition, we exfended EQC fo Gauss-Bonnel gravify oblaining a holonomy correc- 
fed Friedmann equafion, which confains all fhe dynamical informalion of fhe sysfem. Then after 
explaining fhe reconsfrucfion mefhod in holonomy correcfed Gauss-Bonnef gravify, we applied 
our formalism fo cerlain cosmological scenarios, focusing in fhe realization of fhese cosmologies 
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in the context of LQC-corrected F{G) gravity. The cosmological scenarios on which we emphas¬ 
ized are the bouncing cosmologies and the reason for that is that bouncing cosmologies provide 
a quite elegant alternative scenario to the inflationary paradigm. As we evinced, the resulting 
LQC-corrected F{Q) dynamical equations are quite complicated, so we performed a perturbative 
expansion, using as a perturbation parameter, the parameter s which is related to a very important 
physical quantity, the critical density pc, as e = l/pc- Practically, the parameter pc measures how 
quantum is the theory, and the smaller it is, the theory “stretches” in the more quantum era. In the 
perturbation series we used, we assumed that e —)■ 0, so the leading order corrections we found, 
practically quantify the way that the LQC-corrected theory deviates from the classical F{Q) the¬ 
ory. So iteratively, one can recover higher and higher corrections, being though less significant, 
depending on how fast pc tends to infinity. 

Moreover, the possibility of having finite time singularities ll^ and specifically mild types of 
singularities BOll is also quite interesting and should be extensively discussed. 

Finally, the viability of the theories should also be checked and confronted with current ob¬ 
servational data, but this task extends the purpose of this paper, which was the presentation of 
the method of LQC-corrected F{Q) gravity, and its usefulness towards realizing cosmological 
scenarios. We hope however to discuss these topics in the near future. 
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Appendix 


In this appendix, we give the exact form of the parameter C appearing in Eq. (|4^ . It’s detailed 
form is, 


C = 


ClG^ 


+ 


+ 




+ 


+ 


C^G 


+ 


4718592a8 4718592a8 221184a6 147456a6 18432a4 i^2a‘^ 


14155776v/3a9 1572864V3a9 98304v^a7 


294912^30 


CIVG^ Cae^n 


4\/3\/^ 


5Cag3 In 


7 2048V3a5 

4\/3\/^ 


1536v/3a5 l28V3a3 3932160^ 18432a'5 


UC 2 G^ In 


4\/3\/^ 


Ca^ln 


4\/3\/^ C'n^’^/^v/^ln 4\/3\/^ 


1536a3 


16a 


+ 


786432V3, 


a“ 


^9/2y^^ln 4\/3\/^ v/SC'a^^'^^v/^ln 4\/3\/^ 


+ 


+ 


+ 


262144v/3a8 16384a6 

13^^/^\/^ln 4\/3\/^ 3\/3C'|^^/^\/^ln 4\/3\/^ 


49152^31 


+ 


a^ 


1024a4 


^5/2y/^ln 3\/3C'|\/0'\/^lii 4\/3\/^ 


+ 


a^ 


G^ In 


+ 


256^3, 

4\/3\/^ 


+ 


64a^ 


+ —ein 
16 


4\/3\/^ 


■1 2 




131072a6 


+ 




■1 2 


13^2 In 


61440“^ 


+ 


Ay/^yf^ 


■1 2 


512a2 


■1 2 


Ca^^/^v/^ln Ay/^yfofl 2>y/^C2G^^^'/^Aa AypAyf^ 9\/3Ca0^/^v/^ln 4\/3\/^ 


(59) 


■1 2 


131072^3, 


a' 


8192a5 


9\/3Ca\/^'v/^lii 4\/3\/^ 


9 


32a 


H- yf^yfG\/l^\vL Ayf^yf^ 

16 


■1 3 


512a3 
^'^(^\/^ln AyfZyf^ 


■1 3 


+ 


65536V3a6 


-1 3 


3\/30®(^\/^ln AypAyf^ 9\/30^/^\/^ln 4\/3\/^ 
+-^ +- 


■1 3 


4096a4 256a2 

and consequently, in the small G limit, it can be approximated by Eq. (SB- 
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